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ABSTRACT 



Aims. We demonstrate the application of our 3D radiative transfer framework in the model atmosphere code PHOENIX for a number 
of spectrum synthesis calculations for very different conditions. 

Methods. The 3DRT framework discussed in the previous papers of this series was added to our general-purpose model atmosphere 
code PHOENIX/ ID and an extended 3D version PH0ENIX/3D was created. The PH0ENIX/3D code is parallelized via the MPI library 
using a hierarchical domain decomposition and displays very good strong scaling. 

Results. We present the results of several test cases for widely different atmosphere conditions and compare the 3D calculations with 
equivalent ID models to assess the internal accuracy of the 3D modeling. In addition, we show the results for a number of parameter- 
ized 3D structures. 

Conclusions. With presently available computational resources it is possible to solve the full 3D radiative transfer (including scatter- 
ing) problem with the same micro-physics as included in ID modeling. 

Key words. Radiative transfer - Scattering 



1. Introduction 



In 



series 



of 



papers 
(2007); 



lHauschildt & Baronl d2006l); 
lHauschildt & Baronl d2008l 
herj J2009, 



Baron & Hauschildt ._ 

2009h : iBaron. Hauschildt. & Chenl d20()9L hereafter: Papers 
I-V), we have described a framework for the solution of the 
radiative transfer equation in 3D systems (3DRT), including 
a detailed treatment of scattering in continua and lines with a 
non-local operator splitting method. These papers deal solely 
with the radiation transport problem and its numerical solution 
for test cases designed to stress-test the algorithms and codes. 
It is important, however, to apply the radiative transfer codes 
to 'real' problems, e.g., model atmosphere simulations and 
to compare the results to ID equivalents. We have extended 
our general purpose model atmosphere code PHOENIX to use 
the 3DRT framework so that the new version of PHOENIX can 
calculate both ID (PHOENIX/ ID) and 3D (PH0ENIX/3D) models 
and spectra. In this paper we will describe the implementation 
and the results of PHOENIX calculations comparing the results of 
ID and 3D spectrum syntheses for different model parameters. 



2. Method 

In the following discussion we use notation of Papers I-V. The 
basic framework and the methods used for the formal solution 
and the solution of the scattering problem via non-local operator 
splitting are discussed in detail in these papers and will not be 
repeated here. 



3. PH0ENIX/3D implementation and micro-physics 

We have implemented PH0ENIX/3D to use as much as possi- 
ble of the micro-physics of PHOENIX/ ID. This applies to the 
ACES equation of state (Barman, in preparation), to the b-f 
and f-f opacities, to dust opacities, and to the line opacities 
(PH0ENIX/3D is presently restricted to LTE population densi- 
ties). This includes individual line profiles (Gauss profiles for 
weak lines and Voigt profiles for strong lines depending on user- 
selectable selection criteria) for atomic and molecular lines with 
the same physics that is implemented in PHOENIX/ ID, so that the 
results of the opacity calculations are equal for the same physical 
conditions for the two modes of PHOENIX. 

The important considerations of PHOENIX/ 3D implemen- 
tation are memory and CPU time consumption. The memory 
requirements of PH0ENIX/3D compared to PH0ENIX/1D are 
mostly due the the much larger number of voxels in the 3D 
case (typically 10 6 voxels) compared to the ID case (usually 
64-128 layers). As the memory required to store (and to com- 
pute) physical data such as the partial pressures of close to 900 
species or the opacities scales linearly with the number of cells 
(or layers in ID), it is obvious that only very small tests can 
be run without using domain decomposition methods on large 
scale parallel supercomputers. The domain decomposition im- 
plementation of PHOENIX/ 3D distributes the task of solving (and 
storing) the equation of state data and the wavelength depen- 
dent opacities to sets of processes each with its private mem- 
ory. This linearly (with number of processes used) reduces the 
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amount of memory and time required for these tasks. For 1024 
processes, this reduces the memory requirements to just a few 
MB per process to store the full equation of state results. The 
3DRT requires, in comparison, a total of about 450MB for the 
same problem (due to the storage requirements of the non-local 
A* -operator). Including the storage required for the computa- 
tion of the line opacities, this is still just about 0.5GB/process, 
which is small compared to the typically available 4-16GB/core 
(CPU) on modern parallel supercomputers. In order to fully uti- 
lize the available memory per core and to increase flexibility 
we have implemented a hierarchical scheme si milar to the par 



allel PHOENIX/ ID implementation discussed in Hauschildt et al 



1997b and in lBaron & Hauschildt! d 1998b and to the 3DRT paral 



lelization in Paper II: We use a number of 'clusters' of processes 
where every cluster works on a different wavelength. Each clus- 
ter internally uses (on its subset of processes) the domain de- 
composition discussed above and the 3DRT parallelizations dis- 
cussed in Paper I. This scheme can be adjusted to (a) fit the prob- 
lem in the memory available for each core and (b) to optimize 
overall performances (e.g., depending on the number of solid an- 
gle points for the 3DRT solution or the coordinate system used). 
In the calculations presented here, we typically use clusters with 
256-1024 processes, the number of clusters is limited only by 
the number of available CPUs. 

4. Results 

We have calculated a number of test models to compare the 
results of PHOENIX/ID calculations with PH0ENIX/3D results. 
This comparison can be used to adjust the parameters of the 3D 
calculations (number of voxels or solid angle points) to give an 
accuracy that is acceptable for a given investment in computer 
time. The models that we show here were taken from the latest 
PHOENIX/ ID grid (in preparation) of model atmospheres. In all 
stellar models ( ID and 3D) we have used the set of solar abun- 
dances given in lAsplund. Grevesse. & Sau vail (12005b . 



for the different models. In these cases, the error due to the num- 
ber of solid angle points is about 3% and in all tests run the dif- 
ferences between the PHOENIX/ ID fluxes and the F r component 
of the PHOENIX/ 3D calculation is of the same order. The differ- 
ences between the ID and 3D calculations are within the accu- 
racy set by the number of solid angles in the 3D model. In order 
to verify that the errors get smaller with larger number of solid 
angles (as shown in Paper IV for simple test cases), we have run 
test models with 256 2 angles. Three example plots are shown in 
Figs.l8lto[T0lThe results show clearly that the higher solid angle 
resolution reduces the errors in Fg c and F# c considerably and also 
improves the comparison for F r to the ID result, as the higher 
internal accuracy due to more solid angle points also increases 
the internal accuracy of F r , This also shows that in 3D radiative 
transfer calculations the spatial resolution is not the only factor 
governing the quality of the solution, the solid angle resolution 
may in fact be more important, depending on the coordinate sys- 
tem used and the details of the problem that isve calculated a 
number of test models to compare the results of PHOENIX/ ID 
calculations with PH0ENIX/3D results. This comparison can be 
used to adjust the parameters of the 3D calculations (number 
of voxels or solid angle points) to give an accuracy that is ac- 
ceptable for a given investment in computer time. The mod- 
els that we show here were taken from the latest PHOENIX/ ID 
grid (in preparation) of model atmospheres. In all stellar models 
(I D and 3D) we have used the set of solar abundances given 



Asplu nd. Grevesse. & Sauvaj ( 2005 ). We have calculated a 



number of test models to compare the results of PHOENIX/ ID 
calculations with PH0ENIX/3D results. This comparison can be 
used to adjust the parameters of the 3D calculations (number 
of voxels or solid angle points) to give an accuracy that is ac- 
ceptable for a given investment in computer time. The mod- 
els that we show here were taken from the latest PHOENIX/ ID 
grid (in preparation) of model atmospheres. In all stellar models 
(ID and 3D) we have used the set o f solar abundances given in 
lAsplund. Grevesse. & Sauvall (2005). being considered. 



4.1. Stellar Models 

We have computed synthetic spectra for stellar model atmo- 
spheres with the parameters T e ff = 3000K, log(g) = 5.0 (M 
dwarf), r e jf = 5700K, log(g) = 4.5 (solar type star) and 
r eff = 9000K, log(g) = 4.5 (A star). The PHOENIX/ ID mod- 
els were computed with the latest setup in the input physics, 
including the ACES equation of state and the latest version of 
the atomic and molecular line databases. The model structures 
were then used as inputs to PHOENIX/ 3D to calculate synthetic 
spectra with the same sampling rates as the spectra from the 
PHOENIX/ ID calculations. In the PH0ENIX/3D calculations we 
have used a 3D spherical coordinate system with n r = 129, 
ng c = 65 and n$ c = 129 points for a total of about 1M voxels. 
The calculations used (if not specified otherwise) 64 2 solid an- 
gle points. For each object we calculated synthetic spectra with 
PHOENIX/ ID and PHOENIX/ 3D and compare the fluxes of the 
ID spectra to the flux vectors of the 3D results. As in paper IV 
we can use the (0 C , <p c ) components of the 3D flux vector in 3D 
spherical coordinates to estimate the internal accuracy of the so- 
lution (as the Fg c and F^ e components are zero for spherically 
symmetric configurations). Figures [T] to [7] show selected results 



4.2. Scaling 

In order to investigate the strong scaling properties of 
PH0ENIX/3D we have constructed a small test case for a M dwarf 
model with 1000 wavelength points in a 3D spherical coordi- 
nate system with n r = 129, ng c = 65 and n$ c = 129 points 
for a total of about 1M voxels and 64 2 solid angle points and 
ran the calculations with different configurations of the domain 
decomposition and different total numbers of processes. The to- 
tal workload remains constant in these calculations, so this is 
a strong scaling test where the workload per CPU drops as the 
number of processes increase (in contrast to a weak scaling test 
where the workload per process remains constant). The results 
are given in Table Q] In this table, 'n(MPI)' is the total num- 
ber of MPI processes used, 'cluster size' is the number of pro- 
cesses that collaboratively work on a single wavelength (spatial 
domain decomposition) and 'n(cluster)' is the number of such 
clusters, each working on a different wavelength (energy do- 
main decomposition). The product 'cluster size' x 'n(cluster)' 
is always equal to 'n(MPI)'. The column 'Comm' gives the time 
spent in MPI communication to collect the opacities from the 
different processes before the 3DRT calculation starts. The com- 
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munication requirements of the 3DRT calculations are included 
in the 3DRT column. The columns 'line opacity' give the time 
in seconds and scaling efficiency for all line opacity calculations, 
respectively, The columns 'total' give the total time and scaling 
efficiency, respectively, of the overall time spent in the compu- 
tation of the 3D spectrum, this time does not include (small) 
contributions from the EOS solution and the line selection pro- 
cedures. In the largest cluster size of 512 processes each process 
only works on 8 solid angles, whereas in the smallest cluster size 
(128) each process works on 32 solid angles. The work per solid 
angle is not perfectly constant and the amount of communica- 
tion increases linearly as more processes collaborate, therefore, 
the scaling efficiency drops if more than about 512 processes 
are used for this problem size (i.e., number of solid angles). The 
scaling efficiency for the overall problem is quite good, the op- 
timal value is about 98%. The drop-off for cluster sizes of 512 
(and more) is due to (a) the relatively small number of solid an- 
gles leading to very little work for each 3DRT process and rel- 
atively more internal communication time in the 3DRT and (b) 
the small effect of the communication related to the spatial do- 
main decomposition. We could not test setups with more than 
(the maximum available) 2048 processes; however, the test case 
should scale to 256k processes (number of wavelength points 
times cluster size), although for such a setup the overheads for, 
e.g., the solution of the equation of state and the line selection 
would be very noticeable. 

4.3. 3D hydro model of solar convection 

For a test with a computed 3D structure, w e use the same ex- 
ample snapshot structur e from H-G. Ludwig (Caff au et al . 2007; 
Wedemever et al.l 120041) of a radiation-hydrodynamical simula- 
tion of convection in the solar atmosphere as in Paper III. The 
radiation transport calculations were performed with a total of 
141 X 141 x 151 Cartesian grid points in x, y, and z, respectively, 
for a total of 3 002 03 1 voxels, the periodic boundary conditions 
are set in the (horizontal) x, y plane. The 3D radiative transport 
equation is solved for rig = 64 and = 64 solid angle points, so 
that a total of about 12 x 10 9 intensities are calculated for each 
3DRT iteration and wavelength point. For the tests described 
here, we are only using the temperature-pressure structure of 
the hydro model and ignore the velocity field. 

We show example results in Figs.[TT]to[T4]in terms of the x, 
y, and z components of the flux vectors of each outer bound- 
ary voxel. The F z components are, in addition, compared to 
the ID model for the G2V star with the parameters r e ff = 
5700K, log(g) = 4.5 (* symbols in the figures). The general 
shape of the 3D spectra compare well to the ID solar type model, 
of course there are large variations across the horizontal plane. 
In the U V the differences are largest, a number of voxel flux vec- 
tors show strong line emission, whereas the radiative +convective 
equilibrium ID model only shows absorption features. This is to 
be expected as the 3D simulation of convection gives significant 
temperature variations across the volume considered, in partic- 
ular in the horizontal plane. These variations have considerable 
effect on the radiative transfer solution: The horizontal compo- 
nents of the flux vectors of each voxel compared to the length 
of the flux vector F = \F\, F x /F and F y /F, show quite sub- 
stantial variations for different wavelengths. The variations are 



much larger for smaller wavelengths (e.g., in the UV), due to 
the larger temperature dependence of the source functions for 
smaller wavelengths which translates to larger horizontal flux 
components for small wavelengths compared to longer wave- 
lengths. 

The components of the flux vectors in the x — y plane can be 
larger than the z component, strongly dependent on the wave- 
length and on the location of the voxel. This is illustrated in 
Figs. [15]-|22] which show the flowlines of the x—y compo- 
nents of the flux vector at the surface. The flow distances are 
much larger at optical wavelengths than in the UV due to the 
larger UV opacities. The 'pattern' of the horizontal energy flow 
depends strongly on the wavelength, it is also significantly dif- 
ferent in the cores of strong lines compared to the surrounding 
continuum. The horizontal heat exchange could have in turn no- 
ticeable effects on the gas flow pattern. 

4.4. Supernovae 

The modeling of supernova spectra is a very important appli- 
cation of PH0ENIX/3D modeling as it is expected, and explo- 
sion models show, that supernova explosions are intrinsically 3D 
driven. For the calculations shown here we use the Lagrangian 
frame 3DRT in spherical 3D coordinates as discussed in paper 
V. The test model is a simplified model for a type II supernova 
atmosphere with a maximum expansion speed of about 0.13c. 
The model is a simple uniform composition model with the den- 
sity parameterized asp a r , and a "photospheric velocity" of 
vo = 7600 km s , and a model temperature of r mo( j e i = 17000 K. 
These conditions correspond roughly to those of SN 1999em 
seven days after explosion. In Fig. [23] we show the CMF spec- 
trum of the PHOENIX/ 3D run compared to the corresponding 
PHOENIX/ ID synthetic spectrum. Due to computer time limi- 
tations we could only run a relatively small 3D model with 
n r - 129, rig c - 33 and n$ c = 65 points and 128 2 solid an- 
gle points. The small angular resolution causes the scatter in the 
F r plots and the errors in the Fg c and F^ c components. In gen- 
eral the agreement is acceptable for this test run, for a full scale 
3D SN spectrum the resolution in (6 C , (f> c ) should be increased 
to (65, 129) at least and the angular resolution should be at least 
512 2 (which reduces the bandwidth dramatically, see Paper V). 

5. Summary and Conclusions 

We have described first results we have obtained by incorporat- 
ing the 3D radiative transfer framework we have discussed in 
Papers I-V into our general purpose model atmosphere package 
PHOENIX, thus allowing both ID models (PHOENIX/ ID) and 3D 
models (PH0ENIX/3D) with the same micro-physics. We have 
verified and tested PHOENIX/ 3D by computing a number of test 
spectra for ID conditions and comparing the results to the corre- 
sponding PHOENIX/ ID calculations. The conditions range from 
M dwarfs, solar type stars to A stars and Type II supernovae with 
relativistic expansions speeds. In addition, we have calculated 
spectra for a 3D hydrodynamical simulation of solar atmosphere 
convection. These tests demonstrate the it is now possible to 
calculate realistic spectra for 3D configurations including com- 
plex micro-physics. PH0ENIX/3D can be used to calculate syn- 
thetic spectra for a number of complex 3D atmosphere model, 
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including irradiated stars or planets, novae, and supernovae. We 
are currently working on extensions of the 3D radiative transfer 
framework to arbitrary velocity fields in the Euler (for low ve- 
locities, e.g., in convection simulations or planetary winds) and 
the Lagrangian (for Supernovae, accretion disks and matter flow 
in the vicinity of black holes) frames, which will extend the ap- 
plications of PH0ENIX/3D significantly. 
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Table 1. Strong scaling behavior of a M dwarf model test case for different configurations and total number of processors used. See 
text for details. 

Timing: Scaling: 
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Fig. 1. Comparison between the PHOENIX/ ID optical spectrum and the flux vectors across the outermost voxels for the PH0ENIX/3D 
spectra computed for the M dwarf test model ( r e ff = 3000K, log(g) = 5.0, '*' symbols). In the PH0ENIX/3D calculations we have 
used a 3D spherical coordinate system with n r - 129, ng c = 65 and n^ c = 129 points for a total of about 1M voxels. The calculations 
used 64 2 solid angle points. The top panels show the F r component of all outer voxels in linear and logarithmic scales, respectively. 
The bottom panels show the corresponding runs of Fg/F r and F^/F,-, respectively. The should be identically zero and the deviations 
measure the internal accuracy. See Figs. |8]to[T0]for high-accuracy solutions for comparison. The wavelengths are given in A and 
the fluxes are in cgs units. 
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Fig. 2. Comparison between the PHOENIX/ ID near infrared spectrum and the flux vectors across the outermost voxels for the 
PH0ENIX/3D spectra computed for the M dwarf test model ( r eff = 3000K, log(g) = 5.0, '*' symbols). In the PH0ENIX/3D calcula- 
tions we have used a 3D spherical coordinate system with n r = 129, n$ c = 65 and n^ c = 129 points for a total of about 1M voxels. The 
calculations used 64 2 solid angle points. The top panels show the F r component of all outer voxels in linear and logarithmic scales, 
respectively. The bottom panels show the corresponding runs of Fg/F r and F^/Fr, respectively. The should be identically zero and 
the deviations measure the internal accuracy. See Figs. [8] and [10] for high-accuracy solutions for comparison. The wavelengths are 
given in A and the fluxes are in cgs units. 
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Fig. 3. Comparison between the PHOENIX/ ID infrared spectrum and the flux vectors across the outermost voxels for the 
PH0ENIX/3D spectra computed for the M dwarf test model ( r eff = 3000K, log(g) = 5.0, '*' symbols). In the PH0ENIX/3D cal- 
culations we have used a 3D spherical coordinate system with n r = 129, ng c = 65 and n$ c = 129 points for a total of about 
1M voxels. The calculations used 64 2 solid angle points. The top panels show the F r component of all outer voxels in linear and 
logarithmic scales, respectively. The bottom panels show the corresponding runs of Fg/F r and F<p/F r , respectively. The should be 
identically zero and the deviations measure the internal accuracy. See Figs. [8] and [10] for high-accuracy solutions for comparison. 
The wavelengths are given in A and the fluxes are in cgs units. 
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Fig. 4. Comparison between the PHOENIX/ ID near UV spectrum and the flux vectors across the outermost voxels for the 
PH0ENIX/3D spectra computed for the G2V dwarf test model ( r eff = 5700K, log(g) = 4.5, '*' symbols). In the PH0ENIX/3D 
calculations we have used a 3D spherical coordinate system with n r = 129, ng c = 65 and n$ c = 129 points for a total of about 
1M voxels. The calculations used 64 2 solid angle points. The top panels show the F r component of all outer voxels in linear and 
logarithmic scales, respectively. The bottom panels show the corresponding runs of Fg/F r and F<p/F r , respectively. The should be 
identically zero and the deviations measure the internal accuracy. See Figs. [8] and [10] for high-accuracy solutions for comparison. 
The wavelengths are given in A and the fluxes are in cgs units. 
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Fig. 5. Comparison between the PHOENIX/ ID optical spectrum and the flux vectors across the outermost voxels for the PH0ENIX/3D 
spectra computed for the G2V dwarf test model ( r e jf = 5700K, log(g) = 4.5, '*' symbols). In the PH0ENIX/3D calculations we 
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respectively. The bottom panels show the corresponding runs of Fe/F r and F^/Fr, respectively. The should be identically zero and 
the deviations measure the internal accuracy. See Figs. [8] and [10] for high-accuracy solutions for comparison. The wavelengths are 
given in A and the fluxes are in cgs units. 
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Fig. 6. Comparison between the PHOENIX/ ID optical spectrum and the flux vectors across the outermost voxels for the PH0ENIX/3D 
spectra computed for the A dwarf test model ( T e ff = 9000K, log(,g) = 4.5, '*' symbols). In the PH0ENIX/3D calculations we have 
used a 3D spherical coordinate system with n r = 129, ng c = 65 and n$ e = 129 points for a total of about 1M voxels. The calculations 
used 64 2 solid angle points. The top panels show the F r component of all outer voxels in linear and logarithmic scales, respectively. 
The bottom panels show the corresponding runs of Fg/F r and F$/F r , respectively. The should be identically zero and the deviations 
measure the internal accuracy. See Figs.[8]and[T0]for high-accuracy solutions for comparison. The wavelengths are given in A and 
the fluxes are in cgs units. 
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Fig. 7. Comparison between the PHOENIX/ID UV spectrum and the flux vectors across the outermost voxels for the PH0ENIX/3D 
spectra computed for the A dwarf test model ( T e ff = 9000K, log(,g) = 4.5, '*' symbols). In the PH0ENIX/3D calculations we have 
used a 3D spherical coordinate system with n r = 129, ng c = 65 and n^ e = 129 points for a total of about 1M voxels. The calculations 
used 64 2 solid angle points. The top panels show the F r component of all outer voxels in linear and logarithmic scales, respectively. 
The bottom panels show the corresponding runs of Fg/ F r and F,p/F r , respectively. The should be identically zero and the deviations 
measure the internal accuracy. See Figs.[8]and[T0]for high-accuracy solutions for comparison. The wavelengths are given in A and 
the fluxes are in cgs units. 
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Fig. 8. Comparison between the PHOENIX/ ID optical spectrum and the flux vectors across the outermost voxels for the higher 
resolution PH0ENIX/3D spectra computed for the M dwarf test model ( r eff = 3000K, log(g) = 5.0, '*' symbols). In the PH0ENIX/3D 
calculations we have used a 3D spherical coordinate system with n r = 129, tig c = 65 and n$ c = 129 points for a total of about 1M 
voxels. The calculations used 256 2 solid angle points. The top panels show the F r component of all outer voxels in linear and 
logarithmic scales, respectively. The bottom panels show the corresponding runs of Fg/F r and F^jF r , respectively. The should be 
identically zero and the deviations measure the internal accuracy. The wavelengths are given in A and the fluxes are in cgs units. 



Hauschildt and Baron: 3D radiative transfer framework VI 13 




1000 1500 2000 2500 3000 

Wavelength 





0.10 






----- 

□ 


0.05 






% 
J- 


0.00 






2 


-0.05 








-0.10 







1000 1500 2000 2500 3000 

Wavelength 



0.1 OF ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' E 

^ 0.05 E— -E 

X. 0Q T >, ,1 I HHI i ,il I I. - „ ■ „„ ■■,..,>, — ,.- .,,1, 7 

^ -0.05 — — 

-0.1 o E , , , , , , , , , , , , , , , , , , , E 

1000 1500 2000 2500 3000 

Wavelength 

Fig. 9. Comparison between the PHOENIX/ ID UV spectrum and the flux vectors across the outermost voxels for the higher resolution 
PH0ENIX/3D spectra computed for the G2V dwarf test model ( T eS = 5700K, log(g) = 4.5, '*' symbols). In the PH0ENIX/3D 
calculations we have used a 3D spherical coordinate system with n r = 129, ng c = 65 and n$ c = 129 points for a total of about 
1M voxels. The calculations used 256 2 solid angle points. The top panels show the F r component of all outer voxels in linear and 
logarithmic scales, respectively. The bottom panels show the corresponding runs of Fg/F r and F<p/F r , respectively. The should be 
identically zero and the deviations measure the internal accuracy. The wavelengths are given in A and the fluxes are in cgs units. 
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Fig. 10. Comparison between the optical PHOENIX/ ID spectrum and the flux vectors across the outermost voxels for the higher 
resolution PH0ENIX/3D spectra computed for the A dwarf test model ( r eff = 5700K, log(g) = 4.5, '*' symbols). In the PH0ENIX/3D 
calculations we have used a 3D spherical coordinate system with n r = 129, n$ e = 65 and n^ c = 129 points for a total of about 1M 
voxels. The calculations used 256 2 solid angle points. The top panels show the F r component of all outer voxels in linear and 
logarithmic scales, respectively. The bottom panels show the corresponding runs of Fg/F r and F<p/F r , respectively. The should be 
identically zero and the deviations measure the internal accuracy. The wavelengths are given in A and the fluxes are in cgs units. 
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Fig. 11. Comparison between the flux vectors across the outermost voxels for the PH0ENIX/3D UV spectra computed for the 3D 
hydro structure and the PHOENIX/ ID spectrum for the G2V dwarf test model ( T eff = 5700K, log(g) = 4.5). In the PH0ENIX/3D 
calculations we have used a 3D coordinate system with a total of 141 x 141 x 151 Cartesian grid points in x, y, and z, respectively, 
the periodic boundary conditions are set in the (horizontal) x,y plane. The 3D radiative transport equation is solved for ng = 64 and 
= 64 solid angle points. The top panels show the F z component of all outer voxels in linear and logarithmic scales, respectively, 
compared to the results of the ID comparison model. The bottom panels show the corresponding maxima and minima of F J\F\ and 
F y j\F\, respectively, over all surface voxels for each wavelength. These panels show that in the 3D structure even at the surface a 
substantial horizontal energy flow takes place, see also figures[T5]-[22] See text for details. The wavelengths are given in A and the 
fluxes are in cgs units. 
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Fig. 12. Comparison between the flux vectors across the outermost voxels for the PH0ENIX/3D near UV spectra computed for the 
3D hydro structure and the PHOENIX/ ID spectrum for the G2V dwarf test model ( T eB = 5700K, log(g) = 4.5). In the PH0ENIX/3D 
calculations we have used a 3D coordinate system with a total of 141 x 141 x 151 Cartesian grid points in x, y, and z, respectively, 
the periodic boundary conditions are set in the (horizontal) x,y plane. The 3D radiative transport equation is solved for ng - 64 and 
= 64 solid angle points. The top panels show the F z component of all outer voxels in linear and logarithmic scales, respectively, 
compared to the results of the ID comparison model. The bottom panels show the corresponding maxima and minima of F X /\F\ 
and F y /\F\, respectively, over all surface voxels for each wavelength. Note the difference between this result and that shown in the 
bottom two panels of Fig.[TT] See text for details. The wavelengths are given in A and the fluxes are in cgs units. 



Hauschildt and Baron: 3D radiative transfer framework VI 



17 




3500 



3600 



3700 



3800 



3900 



4000 



Wavelength 




3500 




3600 



3700 3800 
Wavelength 



3900 



4000 




3500 



3600 



3700 



3800 



3900 



4000 



Wavelength 




Wavelength 

Fig. 13. Comparison between the flux vectors across the outermost voxels for the PHOENIX/ 3D UV spectra computed for the 3D 
hydro structure and the PHOENIX/ ID spectrum for the G2V dwarf test model ( T eff = 5700K, log(g) = 4.5). In the PH0ENIX/3D 
calculations we have used a 3D coordinate system with a total of 141 x 141 x 151 Cartesian grid points in x, y, and z, respectively, 
the periodic boundary conditions are set in the (horizontal) x,y plane. The 3D radiative transport equation is solved for ng - 64 and 
= 64 solid angle points. The top panels show the F z component of all outer voxels in linear and logarithmic scales, respectively, 
compared to the results of the ID comparison model. The bottom panels show the corresponding maxima and minima of F X /\F\ 
and F y /\F\, respectively, over all surface voxels for each wavelength. Note the difference between this result and that shown in the 
bottom two panels of Fig. QT| See text for details. The wavelengths are given in A and the fluxes are in cgs units. 
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Fig. 14. Comparison between the flux vectors across the outermost voxels for the PHOENIX/ 3D optical spectra computed for the 
3D hydro structure and the PHOENIX/ ID spectrum for the G2V dwarf test model ( T eB = 5700K, log(g) = 4.5). In the PH0ENIX/3D 
calculations we have used a 3D coordinate system with a total of 141 x 141 x 151 Cartesian grid points in x, y, and z, respectively, 
the periodic boundary conditions are set in the (horizontal) x,y plane. The 3D radiative transport equation is solved for ng - 
64 and = 64 solid angle points. The top panels show the F z component of all outer voxels in linear and logarithmic scales, 
respectively, compared to the results of the ID comparison model. The bottom panels show the corresponding runs of F X /\F\ and 
F y /\F\, respectively. See text for details. The wavelengths are given in A and the fluxes are in cgs units. 
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Fig. 15. Illustration of horizontal energy flow for the outermost voxels of the the 3D hydro structure for the red spectral range. The 
graphics shows the flowlines of the x and y components of the flux vector F. Here, a flowline connects points of constant \(F X , F y )\ 
following the direction of (F x , F y ). The 3D radiative transport equation is solved for n g = 64 and = 64 solid angle points. The 
wavelengths are given in A. The normalized x and y voxel coordinates are shown on the x and y axes, respectively. 
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Fig. 16. Illustration of horizontal energy flow for the outermost voxels of the the 3D hydro structure for the red spectral range. The 
graphics shows the flowlines of the x and y components of the flux vector F. Here, a flowline connects points of constant \(F X , F y )\ 
following the direction of (F x , F y ). The 3D radiative transport equation is solved for n g = 64 and = 64 solid angle points. The 
wavelengths are given in A. The normalized x and y voxel coordinates are shown on the x and y axes, respectively. 
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Fig. 17. Illustration of horizontal energy flow for the outermost voxels of the the 3D hydro structure for the red spectral range. The 
graphics shows the flowlines of the x and y components of the flux vector F. Here, a flowline connects points of constant \(F X , F y )\ 
following the direction of (F x , F y ). The 3D radiative transport equation is solved for n g = 64 and = 64 solid angle points. The 
wavelengths are given in A. The normalized x and y voxel coordinates are shown on the x and y axes, respectively. 
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Fig. 18. Illustration of horizontal energy flow for the outermost voxels of the the 3D hydro structure for the red spectral range. The 
graphics shows the flowlines of the x and y components of the flux vector F. Here, a flowline connects points of constant \(F X , F y )\ 
following the direction of (F x , F y ). The 3D radiative transport equation is solved for n g = 64 and = 64 solid angle points. The 
wavelengths are given in A. The normalized x and y voxel coordinates are shown on the x and y axes, respectively. 
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Fig. 19. Illustration of horizontal energy flow for the outermost voxels of the the 3D hydro structure for the red spectral range. The 
graphics shows the flowlines of the x and y components of the flux vector F. Here, a flowline connects points of constant \{F X , F y )\ 
following the direction of (F x , F y ). The 3D radiative transport equation is solved for n g = 64 and = 64 solid angle points. The 
wavelengths are given in A. The normalized x and y voxel coordinates are shown on the x and y axes, respectively. 
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Fig. 20. Illustration of horizontal energy flow for the outermost voxels of the the 3D hydro structure for the red spectral range. The 
graphics shows the flowlines of the x and y components of the flux vector F. Here, a flowline connects points of constant \(F X , F y )\ 
following the direction of (F x , F y ). The 3D radiative transport equation is solved for n g = 64 and = 64 solid angle points. The 
wavelengths are given in A. The normalized x and y voxel coordinates are shown on the x and y axes, respectively. 
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Fig. 21. Illustration of horizontal energy flow for the outermost voxels of the the 3D hydro structure for the red spectral range. The 
graphics shows the flowlines of the x and y components of the flux vector F. Here, a flowline connects points of constant \(F X , F y )\ 
following the direction of (F x , F y ). The 3D radiative transport equation is solved for n g = 64 and = 64 solid angle points. The 
wavelengths are given in A. The normalized x and y voxel coordinates are shown on the x and y axes, respectively. 
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Fig. 22. Illustration of horizontal energy flow for the outermost voxels of the the 3D hydro structure for the red spectral range. The 
graphics shows the flowlines of the x and y components of the flux vector F. Here, a flowline connects points of constant \(F X , F y )\ 
following the direction of (F x , F y ). The 3D radiative transport equation is solved for n g = 64 and = 64 solid angle points. The 
wavelengths are given in A. The normalized x and y voxel coordinates are shown on the x and y axes, respectively. 
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Fig. 23. Comparison between the PHOENIX/ ID co-moving frame UV spectrum (computed with 256 layers, * symbols) and the 
co-moving frame flux vectors across the outermost voxels for the PHOENIX/ 3D spectra computed for the Supernova test model. In 
the PH0ENIX/3D calculations we have used a 3D spherical coordinate system with n r = 129, ng c - 33 and n^, c = 65 points for a total 
of about 275k voxels. The calculations used 128 2 solid angle points. The top panels show the F r component of all outer voxels in 
linear and logarithmic scales, respectively. The bottom panels show the corresponding runs of Fg/F r and F^/Fr, respectively. The 
should be identically zero and the deviations measure the internal accuracy. The wavelengths are given in A and the fluxes are in 
cgs units. 
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